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Spatial structure of X-ray filaments in SN 1006 
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ABSTRACT 

The theory of Non-Linear Diffusive Shock Acceleration (NLDSA) predicts the for- 
mation of a precursor upstream of the shock, where accelerated particles diffuse and 
induce magnetic field amplification through streaming instability. The non detection 
of this precursor in X-rays in Chandra observations of the north-eastern region of SN 
1006 (G329.6+14.6) led to impose an upper limit to the X-ray emission generated by 
accelerated electrons diffusing in this precurs or, at an emissivit y level of < 1.5 per cent 
of the emission from the downstream region ()Long et al. 2003). This has been used as 
an argument against Fermi acceleration at this shock. Here we calculate the spectrum 
and spatial distribution of accelerated particles in SN 1006 and show that Chandra 
results (including more recent data) are in perfect agreement with the predictions 
of NLDSA suggesting efficient particle acceleration and magnetic field amplification 
upstream of the shock by a factor ~ 10. 

Key words: X-ray - supernova remnant - particle acceleration - shock waves - 
magnetic field amplification 
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The blast waves produced by supernova remnants (SNRs) 
expanding in the interstellar medium (ISM) are thought to 
be the sites where the bulk of Galactic cosmic rays (CR) 
are accelerated. The simplest conversion mechanism of ki- 
netic energy of the expanding plasma into accelerated par- 
ticles is Diffusive Shock Acceleration (DSA). Narrow rims 
of non-thermal X-ray emission, associ ated with the outer 
shocks of several young SNRs (see e.g. IParizot et afll2006t 
and references therein), are commonly interpreted as due to 
synchrotron emissio n radiated by ~ TeV electrons (for the 
case of SN 1006 see iBamba et al~ll2003t ). These X-ray fila- 
ments might also represent an indirect evidence of efficient 
acceleration of ions at the same shocks, since their thick- 
ness is most easily explained in terms of severe synchrotron 
losses of elec trons in inten se magnetic fields (in the range 
100-500 juG, iBallet ||200§ ). to be compared with the ISM 
magnetic field, ~ 3/^G. In order to account for this differ- 
ence, it has been proposed that magnetic fields are amplified 
in the shock region by the streaming of cosmic ray ions. This 
gives rise to an interesting scenario in which magnetic fields 
are generated by cosmic rays while in turn making CR ac- 
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celeration more efficient in a complex chain of non linear 
effects. If indeed the magnetic field amplification (MFA) is 
due to CR stre aming, this may happen either via resonan t 
jSkilling I [l975h or non-resonant instabilities (|Bell 1 120041 1 . 
though many variations on the theme are possible. The pe- 
culiarity of CR induced MFA is that it occurs upstream of 
the shock thereby reducing the acceleration time. On the 
other hand MFA could also be generated downs tream by 
fluid instabilities (e.g. iGiacalone fc Jokipiil 120071 ) induced 
by the shock corrugation while propagating in a inhomoge- 
neous medium. Finally the narrow filaments might be due 
to damping of the magnetic field (|Pohl et al.ll200j ) rather 
than to severe synchrotron losses, although correspondingly 
thin fil aments in the radio emi ssion do not seem to be ob- 
served |Rothenflug et al1l2004h . The only way to discrim- 
inate among these different scenarios is through extensive 
comparison between theory and data. 

The spatial profile of the X-ray brightness around the 
shock provides invaluable information on the acceleration 
mechanism and the magnetic field generation. Based on the 
non detection of X-rays with Chandra fro m the upstream re - 
gion of the north-eastern limb of SN 1006. lLong et all l|2003h 
concluded that any X-ray halo upstream must be at least a 
factor ~ 70 fainter than the peak brightness of the shell. 
Based on this finding, these authors suggested that Fermi 
acceleration is not the acceleration mechanism at work or 
that the jump in magnetic field at the shock is 3> 4. They 
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also question the possibility that observations may be ex- 
plained by the non linear theory of DSA. 

In this Letter we apply the theory of NLDSA to show 
that its predictions provide an excellent description of Chan- 
dra observations, which in turn may be used to constrain 
the acceleration efficiency in SN 1006. A p revious attempt 
at us ing the same technique was made in (Bcrczhko et al. 
2003), where a different filament was studied and some as- 
pects of the theory were treated in a more simplified way 
(for instance self-consistent, space dependent MFA in the 
precursor was not accounted for). 



2 NLDSA MODEL FOR SN 1006 

In NLDSA theory, the overall shock structure and the out- 
come of the particle acceleration process are inextricably 
linked. When acceleration is efficient, the pressure of accel- 
erated particles affects the shock dynamics, leading to the 
formation of a precursor, namely a region where the fluid ve- 
locity progressively decreases while approaching the shock 
from far upstream. This causes the total compression fac- 
tor between upstream infinity and downstream to exceed 4. 
At the same time the streaming of accelerated particles is 
responsible for instabilities that lead to MFA. In turn, the 
fluid profile in the precursor and the amplified, turbulent 
magnetic field, with the scattering it provides (and possibly 
the induced energy losses), determine the efficiency of par- 
ticle acceleration and the resulting spectrum, including its 
high energy cutoff. 

Here the shock structure and the accelerate d part i- 
cle spectrum are computed as in iMorlino et all (|20Qgh : 
the basic struct u re of t he ca lculation is that proposed by 
lAmato fc Blasi I (|2005l . | 2006), but the conservation equa- 
tions are modified so as to also take into account the 
dynamical reaction o f the self-generated magnetic field 
l|Caprioli et al. 1 12009^ 1. The compressibility of the overall 
fluid is deeply affected by the magnetic contribution as soon 
as the magnetic field pressure becomes comparable to the 
thermal pressure upstream. The net result is that the total 
compression ratio between downstream and upstream infin- 
ity, Rtot = U0/U2, decreases with respect to the case in which 
this dynamical backreaction is neglected, thereby leading to 
a smoother precursor and accelerated particle spectra that 
are less concave and closer to power-laws. 

The normalization of the proton spectrum is an out- 
put of our non linear calculation, once a recipe for injection 
in the acceleration process is e stablished. Here we follow 
iBlasi. Gabici fc Vannonil l|2005l ). and assume that the frac- 



tion of injected particles, 77i n j, is 



7)inj = 4/ (30F) (R Bub - 1) e" e 



(1) 



where R 3U b = u\/u 2 is the compression ratio at the subshock 
and £ ~ 2 — 4 is the ratio between the injection momentum 
and the momentum of thermal particles downstream, pth,2- 
While £ parametrizes the poorly known microphysics of the 
injection process, Pth, 2 is an output of the calculation: as a 
result, the injection efficiency is affected by the dynamical 
reaction of both particles and fields. 

Energy losses are not important for accelerated protons, 
so that their maximum energy is obtained by equating the 
acceleration time and the age of the remnant, t acc (p Pt max ) = 



tsNR (|Blasi et al.ll2007l) . For electrons the maximum energy 
is reached when the acceleration time equals the minimum 
between the age of the remnant and the loss time scale, 
ti, weighed by the residence times, t r , upstream and down- 
stream of the shock (subscripts 1 and 2 respectively), over 
one cycle: 



tacc(p) = [U,l(p) +tr i2 (p)] 



+ 



Tl tl (B 1 ,p) n i2 (B 2 ,p) 



•(2) 



Once i r ,i and t r>2 are written explicit ly in the contex t 
of NLDSA (from Eqs. (25) and (26) of IBlasi et all 120071 '). 
Eq. ((2| must be solved numerically for p e ,max- The following 
approximate analytical solut ion can be found whe n only syn- 
chrotron losses are relevant (|Morlino et al.| [2009. Eq. (4)): 



= 137 H(p e 



"° 1 TeV ( 3 ) 
3000 km/s J c ' v; 



Here H(p) is a function of the shock modification tak- 
ing into account the mean plasma speed, u p = uoU p (p), 
that a particle with momentum p exp e rience s in the pre- 
c ursor (see Eq. ( 8) of lAmato fe Blasi I (120051 ) and Eq. (4) 
of IMorlino etal] |2009l )). Here H(p) is normalized in such a 
way that H(p) — 1 in the test-particle limit. 

The electron spectrum at the shock, f e .o{p), follows the 
proton spectrum for p <C p e ,max, provided the diffusion co- 
efficient is the same for both species. The normalization of 
the electron spectrum relative to the proton spectrum, K ep , 
is instead unconstrained from the dynamics, since electrons 
have no dynamical role, and can only be derived from ob- 
servations. 

The electron spectrum at energies around and above 
p e ,max, namely the shape of the cutoff, has never been com- 
puted in the context of NLDSA. Since the spectra we find 
for electrons at p < p e , m ax are not far from being power laws 
wi th slope ~ 4, we adopt the modi fication factor calculated 
bv lZirakashvili fc Aharonianl |2007l ) for strong shocks in the 
test particle regime. The resulting electron spectrum at the 
shock, in the loss dominated case, is: 



/ e ,o(p) = K ep f Pi0 (p) 1 + 0.523 (p/p. 



■(4) 



An important point to notice in Eq. @ is that the cut- 
off is super-exponential, which reflects in the shape of the 
synchrotron spectrum, making it different from what is as- 
sumed in most of the previous work on the subject. On the 
other hand, if the maximum momentum is determined by 
the age of the remnant, then the cutoff shape is a simple 
exponential. 

While a few different mechanisms have been proposed 
for MFA at SNR sh ocks, here we f ocus on the waves that are 
resonantly excited (|Skilling II 19751 ) upstream by the stream- 
ing of cosmic rays accelerated at the shock. While stream- 
ing cos mic rays can excite very effectively also non-resonant 
modes ((BeiDEooi), these are likely dominant only at early 
times, while after the beginni ng of the Sedov phase MFA is 
mostly due to resonant waves (|Amato fc Blasi Il2009l ). More- 
over the role of these non resonant modes for scattering the 
particles is yet to be clarified. 

When the predictions of linear theory are extrapolated 
to the non-linear regime of MFA (which one is forced to do 
for lack of a better treatment), the resulting field strengths 
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are in agreement with the values inferred by identifying the 
thickness of the X-ray filaments with the synchrotron loss 
length of the highest energy electrons. The strength of the 
magnetic field at the position x upstream, 8B(x), in the 
absence of damping, can be estimated from the saturation 
cond ition. For modified shocks this reads ijCaprioli et al. I 
120091 . Eq. (43)): 



5B(xf 
87rpo«o 



= U(x) 



-3/2 



U(x) 



4M A ,o 



(5) 



where the Ihs is the turbulent magnetic pressure normal- 
ized to the incoming momentum flux at upstream infinity, 
U(x) = u(x)/uo is the normalized velocity profile in the pre- 
cursor, and Ma,o = Uo/va with va the Alfven velocity at 
upstream infinity, where only the background magnetic field, 
Bo, assumed parallel to the shock normal, enters. Eq. (0 
correctly describes the effect of compression in the shock 
precursor through the term U(x)~ 3 ^ 2 . 

The magnetic field downstream of the subshock is fur- 
ther enhanced by compression, so that B 2 = Rb B\, where 
B\ is the magnetic field immediately upstream of the sub- 
shock and, for Alfven waves (turbulence perpendicular to 
the shock normal), Rb = Rsub- We neglect any kind of mag- 
netic damping both in the upstream and in the downstream 
region. 



3 THE SIZE OF THE X-RAY FILAMENTS 

SN 1006 is one of the best suited remnants for the study of 
the spatial profile of X-ray emission since its non-thermal 
rims, with a typical thickness of 30" -40", are well resolved 
by Chandra, whose angular resolution is Acj> ~ 1". The spa- 
tial resolution corresponding to this value of A(j>, Ax = 
3 x 10 16 (d/2kpc) (A<^>/arcsec) cm on a source at a distance 
d, has to be compared with the extension of the region in 
the upstream (downstream) fluid where electrons radiate at 
the observed frequency uq, i.e. A7? 1 ( 2 )(^o)- Before illustrat- 
ing the results of the detailed calculation, we provide some 
simple but useful approximate relations. 

In the upstream region electrons diffuse up to a dis- 
tance Ldiff = D(p) /u p (p) from the shock, where D is the 
diffusion coefficient and u p is the effective fluid velocity ex- 
perienced by particles with momentum p. The spatial dis- 
tribution of electrons, /i(p), is roughly c onstant up t o the 
diffusion length Ldiff, an d is then cut off (jBlasi 1120091 ). As- 
suming Bohm diffusion, Db = p<? / (3eB) , we find: 



ARi(v) ~ Ldift (H 



D B (p) 



pc 



x 10 1 



IkeV/ 



u p (p) 3eBi Up (p) 
Bi V* f uoU P {p) 



10 uG 



3000 km/s 



(6) 



where we used the fact that synchrotron emission of elec- 
trons with energy E in a fiG field peaks at energy hv = 
1.5 x 10 _2 L? m g i?TeV 2 eV. Since Ldiff strongly depends on the 
magnetic field strength, the possibility of resolving the up- 
stream X-ray halo depends on the level of upstream MFA. 

In the upstream region, electrons with E — B max radi- 
ate at a typical frequency 



where we used the expression in Eq. ([3]) for L max . Clearly 
Eq. ((6J holds only for E < E m ax , since at higher energies the 
diffusion region upstream becomes regulated by synchrotron 
losses: 



ARi(v >ui)~ Ldifi (hu!) 



(8) 



The size of the radiating region downstream is determined 
by the competition between advection and diffusion of the 
electrons during their life-time, namely L^ = ui t\ (p) oc 
p _1 and L^ifl = \J D(p) n(p) (the latter is independent 
of energy for Bohm diffusion). These two length-scales are 
equal for particle momenta corresponding to photon energy: 



f up y fRtot\ 

UOOOkms- 1 / V 4 J 



keV 



hv 2 = 1.00 , 

\ 3000 km s- 

At frequencies lower than v 2 : 



\ r, / , \ r loss u 1 \ —1/2 

AR 2 (y < v 2 ) ~ L adv = — — Ti(p)oci> 

tttot 



(9) 



(10) 



while for v > v 2 the downstream emission comes from a 
region of constant size: 



loss , r w ln 17 / B2 



AR 2 {u > i/ 2 ) ~ LdTff = 4.5 x 10 



V 40pG 



(11) 



In order to compare ARi and AR 2 with the ob- 
served rim thickness, we need to take into account projec- 
tion effects. Let us define the observed thickness, AR°^, 
as the width over which the emissivity drops to half the 
peak value. Assuming spherical geometry, an emission pro- 
file I v oc exp[— x/AR] and ARu^) much smaller than the 
SNR radius, the observed widths are ARf s = 0.7ABi 
and Ai?2 bs = 4.6AJ?2- Hence projection effects make the 
observed width of the downstream emission region appear 
larger than the physical one. 

When we observe the filaments at frequencies v > v\ % v% % 
the ratio between the observed thickness downstream and 
upsteam can be estimated as: 



AR% 



4-6 L\°{§ 



AR° bs 0.7 L d iff (p e 



12 



1 + RsRtotUp 
RKl-R^Up 1 ) 



(12) 



hv x = 0.28#(p e , max )(u /3000km s" 



keV, 



(7) 



where the last equality is obtained using the full expression 
for H (p) that appears in Eq.[g](see Eq. (4) of iMorlino et all 
2009). Remarkably this result is independent of magnetic 
field strength, shock speed and electron energy, while it only 
depends on the compression ratios (magnetic compression at 
the subshock Rb and total compression Rtot) and on U p = 
Up/u . Since Rsub /Rtot < Up < 1 and \/TT ^ R B < 4, using 
reasonable values for Rtot the ratio A7?2 bs /Aiii bs is always 
between 6 and 8. Since the rim in SN 1006 has an overall 
width of 30" -40" for v ^ Ui, u 2 , the result above implies that 
the upstream halo should have an observed width of < 6". 



4 COMPARISON BETWEEN THEORY AND 
OBSERVATIONS 

Our model of NLDSA contains four free parameters, namely: 
the density of the unshocked ISM, no, the shock speed, 
wo, the proton injection efficiency, £, and the relative nor- 
malization between electron and proton spectra, K ep . SN 
1006 is thought to be the remnant of a type la SN 
explosion and to expand in a fairly uniform ISM. We 
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Figure 1. Chandra X-ray image of the northeastern limb showed 
in squared root color scale. The white rectangle (50" X 120") is the 
region used to extract the X-ray radial profile plotted in Fig. [2] 
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Figure 2. Radial profile of the emission in the 0.8-2.0 keV band, 
extracted from the rectangular region in Fig. [T] (f/iin black line). 
Overplotted are the theoretical predictions of NLDSA for 3 dif- 
ferent injection efficiencies (thick lines) and those of test-particle 
theory for two different values of the pre-existing turbulent mag- 
netic field upstream (thin lines). 
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Bj (jiiG) 


B 2 (jiG) 




Rtot 


3.90 


47 


175 


3.78 


6.39 


3.95 


23 


90 


3.93 


5.53 


4.15 


5.3 


21 


4.00 


4.08 


CO 


3.0 


10 


4.00 


4.00 


oo 


10 


33 


4.00 


4.00 



of the proper motion of X-ray filaments by iKatsuda et al.l 
< |2009h . which give u = 4670 (d/2kpc) km/s. 

We extract the X-ray emission profile from recent Chan- 
dra observations (20 08 June 24; ObsID.9107; PI: Petre; see 
IKatsuda et al.l [20091 ) . Fig. [T] shows the Chandra image of 
the northeastern limb. We selected the rectangular region 
of 50" x 120" highlighted by the white box in the figure, 
so as to exclude regions where multiple shocks overlap. The 
resulting radial profile, in the energy band 0.8-2.0 keV, is 
plotted as a histogram i n Fig . [5] and is in good agreement 
with that of lLong et all i|2003l ). 

The curves in Fig. [2] refer to different values of the in- 
jection parameter £ = 3.90, 3.95, 4.15 and Bo = 3/j,G. The 
two curves labelled as B\ = 3/iG and 10/j.G are test particle 
cases with a turbulent magnetic field as indicated, and with 
the perpendicular components compressed at the shock, so 
that £?2 = vTTBi. In the other cases the magnetic field at 
a position x in the precursor is calculated as described by 
Eq. ([5|, and then compressed at the subshock, resulting in 
the values of B\ and B2 reported in Table HJ together with 
the compression factors R su b and Rtot- 

It is clear that DSA in the test particle regime cannot 
explain either the radial profile in the precursor or that in 
the downstream region, although the latter might still be 
explained by requiring some magnetic damping. The inter- 
mediate case, £ = 4.15, which corresponds to modest in- 
jection, is also inappropriate for the same reasons. More- 
over, the test particle and the inefficient cases show a poten- 
tially detectable precursor, which is however not observed. 
Remarkably, the solution with £ = 3.95 (corresponding to 
B2 — 90/j.G) fits well the downstream X-ray profile and 
shows a modest upstream halo that is however fully con- 
sistent with the observations. For larger injection efficien- 
cies the upstream halo becomes smaller than the Chandra 
resolution, but the predicted downstream rim is too narrow 
compared with observations. 



adopt np = 0.05 cm" 3 , based on the measurement by 
lAcero et al. I (|2007l ) (also consistent with a hadronic inter- 
pre tation of the TeV emission recently det ected by H.E.S.S.; 
see iBerezhko. Ksenofontov fc Volk 1120091 ). We estimate the 
shock speed using a code for the evolution of type la SNRs 
ijCaprioli. Amato fe Blasil I2010T I. fixing the explosion en- 
ergy and the mass of the ejecta to the standard values 
Esn = 10 51 erg and M e j e = 1.4M Q . Once rio is fixed, given 
the known age of the remnant, we estimate no = 4330km/s 
and a radius Rsnr = 7.7 pc. Comparing this latter value 
with the observed angular diameter of 30', the inferred dis- 
tance of SN 1006 is d = 1.77 kpc. This makes our estimate 
of uq in remarkable agreement with recent measurements 



5 DISCUSSION AND CONCLUSIONS 

We studied the observed radial profile of the X-ray emis- 
sion in SN 1006 in the context of the theory of NLDSA. 
We concentrated our attention on a filament of X-ray emis- 
sion which was already studied bv lLong et al.l ((2003) since 
it shows a pronounced jump in emission by a factor ~ 70 
at what appears to be the shock surface. This filament is 
rather clean from the observational point of view, in that 
there seems to be only one shock surface within the observed 
slice. 

The conclusions of our investigation can be summarized 
as follows: 
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1) The radial profile of the X-ray emission is very well 
described within the context of NLDSA, as a consequence of 
a large magnetic field amplification upstream of the shock. 
This scenario describes both the downstream narrow fila- 
ment and the non detection of a precursor upstream of the 
shock. The main reason for this non-detection is that the 
magnetic field upstream becomes strong enough to make 
the diffusion/loss region of electrons comparable with the 
angular resolution of Chandra. Moreover, if magnetic field 
is amplified upstream through the excitation of a stream- 
ing instability, the gradient in the precursor magnetic field 
makes the size of the X-ray emission even smaller. 

Fixing no = 0.05cm -3 , we best fit the observed ra- 
dial profile with £ = 3.95 which implies a CR accelera- 
tion efficiency of ~ 29 per cent and an electron/proton 
ratio K Bp ~ 1.3 10~ 4 . An inefficient scenario of CR accel- 
eration appears to be inadequate to fit observations: the 
narrow filaments would require some level of magnetic field 
damping downstream, which however does not seem to be 
compatible with the ab sence of filaments in the radio band 
|Rothenflug et all 120041 1. Upstream, the inefficient scenario 
typically leads to broader emission regions than observed by 
Chandra. 

We stress that our conclusions do not strongly depend 
on the assumed values of the upstream magnetic field and 
shock speed. An equally good fit of the X-ray radial profile 
can be obtained for Bo — 5uG (about the highest reasonable 
value at the galactic latitude of SN1006) and for a shock 
velocity of 4670 km/s corresponding to a larger distance 
of the remnant (d = 2kpc) and to an exposion energy of 
1.5 • 10 51 erg. For a magnetic field strength Bo = 5uG the 
best-fit injection parameter needs to be changed (£ = 3.98), 
corresponding to a slightly smaller acceleration efficiency. 
For a shock velocity of 4670 km/s a fit can be obtained for 
£ = 3.93. 

2) O ur conclusions are in disagreement with the point 
made by lLong et al.l (|2003h that one needs a jump in the 
magnetic field at the shock which is 2> 4, equivalent to MFA 
downstream. We find that the important in gredient is that 
there is MFA upstream, as also found by iBerezhko et al.l 
|2003t ). although adopting values of the ambient density and 
shock velocity that have now been excluded by observations 
|Acero et al. II2007T ). Our calculations show that in t h e case 
of SN 1006, inclusion of MFA a la iBell fe Lucekl (|200lh 
(which would further enhance the magnetic field upstream) 
is not necessary. 

3) Contrary to IBerezhko et al.l l |2003t ) we think that an 
efficient CR acceleration is not the only reasonably thinkable 
condition to explain the X-ray observations. The non detec- 
tion of a CR precursor upstream can also be attributed, at 
least qualitatively, to the presence of a quasi-perpendicular 
magnetic field, as was also inferred bv lBamba et al.l (|2003l ') 
(even if their conclusions were based on the incorrect as- 
sumption of no magnetic jump at the shock surface); in this 
case CRs can be accelerated to rather high energies l|jokroii| 
1 19871 ). In this situation the magnetic field could be ampli- 
fied downstream of the shock as due to fluid instabilities 
if the shock propagates i n a inhomogeneous medium (e.g. 
iGiacalone fe Jokroiil [2007h . Whether this configuration can 
explain the details of the X-ray radial profile should be sub- 
ject of a dedicated investigation. 
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